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We describe a Josephson device composed of two superconductors separated by two interacting 
quantum dots in parallel, as a probe for Cooper pair splitting. In addition to sequential tunneling of 
electrons through each dot, an additional transport channel exists in this system: crossed Andreev 
reflection, where a Cooper pair from the source is split between the two dots and recombined in 
the drain superconductor. Unlike non-equilibrium scenarios for Cooper pair splitting which involves 
superconducting/normal metal “forks”, our proposal relies on an Aharonov-Bohm measurement of 
the DC Josephson current when a flux is inserted between the two dots. We provide a path integral 
approach to treat arbitrary transparencies, and we explore all contributions for the individual phases 
(0 or 7r) of the quantum dots. We propose a definition of the Cooper pair splitting efficiency for arbi¬ 
trary transparencies, which allows us to find the phase associations which favor the crossed Andreev 
process. Possible applications to experiments using nanowires as quantum dots are discussed. 

PACS numbers: 74.50.-tr 85.25.Dq 73.21.La 


I. INTRODUCTION 

In superconductors, Cooper pairs provide a natural 
source of entanglement 1 in both spin and momentum 
space. If for some reason the constituent electrons of 
this pair are separated and are allowed to propagate in 
two different metallic leads or nanodevices such as quan¬ 
tum dots, one expects that the entanglement is preserved 
because the tunneling processes are spin preserving. This 
process is called cross Andreev reflection (CAR), and has 
been the focus of both theoretical and experimental in¬ 
vestigations in the last two decades. The initial man¬ 
ifestation of CAR was proposed theoretically for non¬ 
local current as well as for non-equilibrium noise 
cross-correlations (current-current fluctuations between 
the superconductor and the two outgoing devices).®^ 
Indeed, the scattering formalism of electron and hole 
transport showed that noise cross-correlations could be 
positive. Strictly speaking the positive cross-correlation 
signal does not constitute a rigorous proof of electrons 
entanglement originating from superconductors, but it 
certainly provides evidence of Cooper pair splitting. Al¬ 
ternative approaches exploiting electron energy filtering 
using Coulomb blockade confirmed at the same time the 
nonlocal spin singlet nature of electron pairs in opposite 
quantum dots placed in the vici nity o f a superconductor, 
exploiting T matrix calculations “23 

The possibility of generating nonlocal entangled pairs 
of electrons in condensed matter settings bears funda¬ 
mental applications in the context of quantum informa¬ 
tion theory. Tests of quantum entanglement followed 
these works, based on the idea that Bell inequality viola¬ 
tion measurements could be implemented via noise cross¬ 
correlations wit h the so-called superconducting Cooper 
pair splitter. 17 ^ Moreover, these ideas have b een a p- 
plied to the paradigm of quantum teleportation] 19 l 2 ° l as 
well as in other applications j 21 l 22 l 

On the experimental side, attention has mainly fo¬ 


cused on nonlocal current measurement on the Cooper 
beam splitterpS a device where typically the supercon¬ 
ducting source of electrons is connected to two leads, 
sometimes via embedded quantum dots.^^ Under spe¬ 
cific gate voltages imposed on such dots, it is possible 
to trigger electronic transport in the two outgoing con¬ 
ductors. Only a single experiment managed to measure 
positive noise cross-correlations when the source of elec¬ 
trons was rendered superconducting. 28 

The main challenge with these proposals resides in the 
fact that they rely on non-equilibrium measurements. 
Strictly speaking, these measurements, although to be 
saluted, constitute only an indirect evidence of Cooper 
pair splitting, while noise cross correlation measurements 
represent a considerable ordeal due to the poor signal to 
noise ratio, and no attempt has been tried so far to re¬ 
produce them. 

A seminal theoretical work has been suggested early on 
to circumvent these difficulties by proposing a Joseph¬ 
son equilibrium current geometry to test Cooper pair 
splitting.^ It describes two superconductors (with ap¬ 
plied phase difference) separated by two quantum dots 
placed in parallel. When a Cooper pair is transmitted 
from one superconductor to the other, the two electrons 
can either pass both through a given dot, or they can 
transit through different dots (cf. Fig. [I]). This indeed 
realizes an Aharonov-Bohm (AB) experiment with su¬ 
perconductors as source and drain, driven by an applied 
phase difference. The critical current as a function of the 
AB flux should be 7r periodic if electrons are not split be¬ 
tween the two dots, and 2ir periodic if Cooper pair split¬ 
ting is effective. The clear originality of this proposal 
resides in the fact that unlike non-equilibrium noise se¬ 
tups, here Cooper pair splitting is uncovered using a cur¬ 
rent measurement at equilibrium albeit in a Josephson 
geometry. The calculation was performed perturbatively 
in the tunneling Hamiltonian, with infinite repulsion on 
the dots. Dot gate voltages insured that on average each 
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Figure 1. Illustration of the three different possibilities for the 
transmission of a Cooper pair from the left superconductor to 
the right one. The two electrons of the pair can either be 
both transmitted through the upper dot (top left), through 
the lower dot (top right), or the Cooper pair can be split with 
one electron transmitted through each dot (bottom). 


dot was occupied by a single electron. A complementary 
study appeared a decade later with the same setup 30 and 
perturbative results for dot levels which were assumed to 
be above the superconducting chemical potentials and 
with a finite Coulomb repulsion were also presented. 

First, so far no analysis of this superconducting 
Aharonov-Bohm effect has allowed to go beyond lowest 
order perturbation theory in the tunneling Hamiltonian. 
Advances in superconducting device fabrication 311 seem 
to indicate that by burying nanowires underneath super¬ 
conductors, large transmission can be achieved between 
the resulting quantum dot and the lead. Treating the 
tunnel coupling to all orders of perturbation theory thus 
constitutes a first motivation of our study. 

Secondly, it is established that when a quantum dot is 
embedded in a Josephson junction, away from the Kondo 
regime, the strength of the on-site repulsion, the coupling 
to the leads and the level position of the dot determine 
whether it constitutes a 0 or a it junction (positive or 
negative Josephson amplitude). This has been analyzed 
perturbativelyP^l as well as with path integral formalism, 
followed by a saddle point approximation.^ This latter 
work allows to distinguish between three phases of the 
Josephson junction: a) the tt phase (where the dot is 
singly occupied); b) the Ch 0 ' phase (where it is unoccu¬ 
pied); c) and the 0^ phase (with double occupation). 
These predictions have been verified experimentally a 
decade agop3 In the Josephson setup with two dots in 
parallel which is studied here, each dot can either be in 
the CK 0 ), the C)( 2 ) or the tt phase. In the works of Ref. £39; 
the two electrons are both in a tt junction configuration 
while the work of Ref. [30] have also considered the case 
where both dots are in the 0*-°) phase. However, there is 
to this date no systematic or comparative study specify¬ 
ing which combination of the dot phases may enhance or 
reduce the AB signal for maximal observation, even less 
so for arbitrarily large transparencies. 

The two above points constitute the main motivations 
of the present paper. In this work, we employ the path 


integral formalism to model the AB setup without any 
restrictions on the transmission properties of the sam¬ 
ple. Indeed, we provide results both in the tunneling and 
the high transparency regimes, and we propose a way to 
measure the efficiency of CAR processes in both cases, 
by analyzing the AB signal. We find that the CAR pro¬ 
cesses are optimized when the two quantum dots are in 
the same phase. 

The paper is organized as follows. In Sec. [TT] we pro¬ 
vide a description of the device and of the model with 
which we describe it, and we give an expression of the 
partition function in terms of Grassmann variables. The 
free energy used to derive the Josephson current of this 
nanoSQUID in a non-perturbative manner is presented 
in Sec. m In Sec. m we discuss the possible occupancy 
states of the dots. We propose the definition of a Cooper 
pair splitting efficiency in Sec. |V| which can be computed 
for arbitrary transparencies of the studied junctions. The 
AB signals are calculated first in the tunneling regime 
in Sec. J3 and then in the non-perturbative regime in 
Sec. |VII and all possible phase associations for the dots 
are considered. We discuss our results in Sec. Em] 


II. MODEL AND PARTITION FUNCTION 

The device is illustrated in Fig. [2] For simplicity, the 
two leads consist of the same superconducting material 
with a controllable phase difference which can either be 
imposed by closing the device in a loop geometry or em¬ 
bedding the device in a macroscopic SQUID!2212SI Two 
quantum dots are placed in parallel in the nanogap be¬ 
tween the two electrodes, and a magnetic flux (which is in 
principle independent from the one imposed to trigger a 
DC Josephson signal between the electrodes) threads the 
area between the two dots. Electrons can tunnel from the 
source electrode to the upper or lower dot, but on site 
Coulomb repulsion favors zero, or single occupation on 
the latter. In the presence of a magnetic flux between the 
two dots, because of the AB effect, different phase shifts 
are expected between the two paths that an electron can 
follow to reach the drain electrode. If the separation be¬ 
tween injection points is larger than the coherence length, 
Cooper pairs as a whole pass either by the upper or lower 
path, as in the work of Ref. [3Tj If the injection points 
are closer together than the superconducting coherence 
length however, the two electrons of the pair can travel 
through opposite dots, realizing a CAR process, which 
gives a third contribution to the Josephson current. By 
adjusting the quantum dot energy levels and Coulomb 
interaction, we expect to filter the electrons and eventu¬ 
ally favor the CAR process. Granted, we cannot claim to 
directly measure the degree of entanglement of this CAR 
contribution (like in a Bell inequalities measurement), 
nevertheless, the recombination of these two electrons on 
the destination superconductor could not be possible if 
this outgoing state did not form a Cooper pair, i.e. an 
entangled state. 
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We denote by dP aa the creation operator for an electron 
with spin a =f, on the quantum dot a = U,D and by 
tf’jkcr the creation operator for an electron with momen¬ 
tum k and spin a =f, i in the superconductor j = L, R. 
It is convenient to introduce the Narnbu spinors 

and (1) 

fj i (i = x, y , z) are the Pauli matrices that act in Nambu 
space. The Hamiltonian of the double Josephson junction 
reads 

n= X H °+ E H J + Hf (2) 

a=U,D j—L,R 

H a is the Hamiltonian of the quantum dot a = U, D, 
characterized by its energy level e a and its on-site 
Coulomb repulsion U a - 

Ha — X aa + U a Tl a jTl a ^ , (3) 

<7=t4 

with n arT = cr aa d a a the dot occupation operator per spin. 
Hj is the Hamiltonian of the superconductor j = L , R, 
with gap A and chemical potential y: 

- k 2 

H j = X, %k (& + A °x) V’jfc! ^ k= 2m~ ^ 

k 

Here A is real, the phase difference between electrodes 
has been gauged out and instead appears in the tunnel¬ 
ing Hamiltonian. If we denote by rj a the location of the 
injection point from lead j to dot a, the tunneling Hamil¬ 
tonian H t reads 

H t = X ° lk rja tyk T ja d a + h.c. (5) 

jka 

The tunneling matrices involved in Eq. § read 

Tlu =t L cr z e +l!L ^ ,Jz , Tru =t R <j z e~ l± ^' 7z , (6a) 

Tld =t L a z e +t!t r Lr7z , Trd = t R a z e ~ l ^ az . (6b) 

4> is the phase difference between the superconductors 
while a is related to the magnetic flux $ inside the 


SQUID loop: a = 2n where $o = h/e is the flux quan¬ 
tum. For clarity, the phase shifts acquired by tunneling 
electrons are indicated in Fig. [2] 

We employ a path integral approach in the Matsubara 
formalism in order to compute the partition function of 
the device. We introduce then the eigenvalues of the 
annihilation operators ip jka and d aa written as ipjka and 
d arT respectively. These are Grassmann variables and we 
consider also their conjugates ipjka and d a a as well as the 
collections in Nambu spinors d a and ipjk defined in the 
same way as Eq. |l]). 

The partition function is given by a functional integra¬ 
tion over paths that are /3 antiperiodic: 

Z= J V [d,d,ip,ip] exp [S E (d,d,ip,ip)] . (7) 

da(/3) = -da(0) 

■0jfc(/3)=-bjfc(O) 

The Euclidean action Sr reads 

0 

S E (d,d,ip,ip) = J dr j'H (d, d, ip, ip) 

o 

+ X dad T d a + x; i’jk dripjk i (8) 

a jk ' 

where the matrix elements of the Hamiltonian can be 
written as 

H (d, d, Ip, Ip) = X H a (d a , d a ) + X H jk (ipjk, ipjk) 

a jk 

+Y. H * ,jka ( d a ,d a ,ip jf,,ipjk )• (9) 

jka 

The expressions of Hjk and H t jka are readily obtained 
from Eqs. (4])-([5]) by substituting the annihilation opera¬ 
tors a by their eigenvalues a and the corresponding cre¬ 
ation operators aP by the conjugate Grassmann variables 
a. For the quantum dots, we can also find an expression 
in terms of Nambu spinors as follows 

- - U - 2 

Ha (da, da) — Ca T &a d a & z d a (d a d a ) , (10) 

with E a = £ a + \- 

III. FREE ENERGY AND JOSEPHSON 
CURRENT 

As the lead degrees of freedom are quadratic in the 
Hamiltonian, they can be easily integrated out. The par¬ 
tition function is then expressed as a functional integral 
over the dot Grassmann variables: 

Z = C\ J W V {da, da] exp [~S eS (d, d)] , (11) 

da 08 )=-d„(O) ° 

where ci is the determinant arising from the integration 
of the lead variables, which is independent of <p and a. 
The effective action in Eq. reads 
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S eS (d, d) = ^2 / dr 


£a + 


d a {r) [d T t 2 + £ a <r z \ d a (r) - -y- (d a (r) d a (r)^ j - ^ j dr J dr' d a (r) S a6 (r - r') <4(r') 

J a ’ b o o 

( 12 ) 


where the self-energy term 

E“ & (r) = J2^ (rib - r ^rlG k (T)T ]b (13) 

jk 


involves the Green function of the leads G k (r) which ver¬ 
ifies 


[<9rl 2 + CfcCT z + Aer x ] G fc (r) = <5(t)1 2 . 


(14) 


The quartic terms (d a d a ) 2 in Eq. (121 prohibit an exact 
computation of the partition function. As in Ref. l33l 
we use a Hubbard-Stratonovich transformation to treat 


these terms and we neglect the temporal fluctuations of 
the two auxiliary fields Xjj and Xjj which are introduced: 



Jdr(d a d a ) 2 



dX a e 


~2§~ x a+ x a J dr d a d a 


(15) 


Because both the Green functions G k and the Nambu 
spinor components d aa are fd antiperiodic, we use a Mat- 
subara series expansion F(t) = e~ lulpT F( uj p ) over 

the frequencies co p = (p + |) 2t: / fd. 


Rather than keeping track of a cumbersome device-specific position dependence of the self-energy, we choose to 
introduce a phenomenological parameter rj. It extrapolates between the two most relevant cases: rj = 0 for infinitely 
distant injection points (much more separated than the superconducting coherence length) and 77 = 1 for coinciding 
injection points (much closer than the superconducting coherence length). The dots are now integrated out and the 
partition function becomes 


Z v {(l),a) = c ic 2 



✓»+oo 


d X £> exp 


C-HS.p 

D eff 


(Xu, X D , <j>, a) 


where c 2 is a fermionic determinant arising from the dot variables integration. The effective action reads 


‘-’eff 


(Xu,XD,<f>,a) = 


1 £q + 


2 U a 


X„ ) — 2 


pen 


ln^ 4 det [M^(Xu,X D ,(f>,a)] ), 


(16) 


(17) 


Mp (Xu,X D ,(j), a) = 


- (iuip + Xu) 1 2 +EU&Z- A p (4> - a) 


~B p (<t>, +«) 

— (iljjp + Xu) I 2 + Ed &Z ~ Ap(0 + a) 


(18) 


Ap(<t>) = 


B2((/>,a) = rj 


iio p 1 2 — A ( cos — cr x + 7 sin — a l 


A 2 +w2 


fa . a \ / (j> . <j> \ 

vjj p (^cos - 1 2 + sin — a z J — A I cos - a x + 7 sm ^ a v ) 


(19) 

( 20 ) 


The decay rate is defined as T = 7n/(0)(f| +f|j) where ^(0) is the density of states of the leads at the Fermi level. 
The contact asymmetry is given by 7 = (t\ — + Ir)- 


The CAR process is due to the off-diagonal terms B p 
of the matrices M. p the determinants of which we have 
to compute as a result of the Gaussian integrals over 
quantum dot degrees of freedom. In practice, these off- 
diagonal terms depend on the separation between injec¬ 


tion points R = |rju — r jD|- They have an exponen¬ 
tial decay on the scale of the superconducting coherence 
length. With the microscopic tunneling Hamiltonian for¬ 
mulation of Eq. ([ 5 ]), they also bear fast oscillations on 
the Fermi wavelength, and possibly power law decay de- 
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pending on the dimensionality [rj(R) = (sinfcFi?)/(fcpi?) 
for 3D superconductor^ 35 anc j n0 p 0wer l aw decay 
for quasi-one dimensional superconductor^^ . Existing 
CAR experiment^®' based on nanowire quantum dots 
embedded in the superconducting leads are typically per¬ 
formed for an injection separation which is less than 
the superconducting coherence length. Yet these experi¬ 
ments find no evidence of either Fermi wavelength oscil¬ 
lations or power law decay at all: the nonlocal signal is 
strong and can be optimized by tuning the dot gate volt¬ 
ages. This can be attributed to the proximity effect from 
the bulk superconductors, which acts on the nanowires 
used in the experiments. To keep our discussion more 
general, and avoid such device-specific complications, in 
this work, we use i] as a phenomenological parameter, as 
in Refs. [301 and 133 Most of the results will be displayed 
for the extreme values 77 = 0 and 77 = 1 , but in order 
to show the evolution of the AB signals, we sometimes 
allow it to vary smoothly between these two values. 


To evaluate the partition function of Eq. (16), 
we use a saddle-point method.^ The effective action 
(Xu,Xr>,(f>,a) is computed numerically by sum¬ 
ming over Matsubara frequencies (up to a cut-off much 
larger than the superconducting gap). Its minimum, lo¬ 
cated in [X^j{(f), a), Xp((j), a)], is obtained with a gradi¬ 
ent descent method for fixed <j> and a (the number of 
starting points of the algorithm depends on the sym¬ 
metry of the function to be minimized). The free 
energy is then defined from this minimum value as 
S^’ v (XZ(</> t a),X* D {<l>,a),<j>,a) = pF^^a). The cur¬ 
rent is finally obtained by differentiating the free energy 
with respect to the phase difference <fi: 


J TI {<t>,a) = 2 0^(0, a). ( 21 ) 


The critical current, function of the flux a, is defined as 


J^(cc) = max^ | a) | . 


( 22 ) 


For 77 = 0, it is 7 T periodic. Indeed, in this case, there 
are no CAR processes and for a magnetic flux a, the 
current characteristic (current as a function of the phase 
difference <j>) through the double junction is simply the 
sum of the current characteristics across the independent 
single junctions, one being shifted by —a, the other one 
by +a. Shifting a by 7 r results in adding a phase shift 
— 7 r for one of the two junctions [A p ((f> — a) —> A p {<j> — 
a — 7 r) in the top left block of Eq. (18)] and a phase 
shift -t- 7 r for the other one [Ap{(f> + a) —» A p ((j> + a + 7 r) 
in the bottom right block of Eq. ©]. Using the 27 t 
periodicity in (f> of the current through a single junction, 
and taking the maximum on <fi shows immediately that 
the critical current in unchanged when a —► a+ 7 r. This is 
in contrast to the case with rj = 1 , i.e. in the presence of 
CAR processes, where the off-diagonal terms of Eq. (18) 


implies that only the 27 t periodicity of the critical current 
can be observed. 



Figure 3. Mean occupation number diagram of the quantum 
dot of a single Josephson junction for symmetric couplings 
(tL = ti?) and for F = A, /3 = 50A” 1 . 


IV. 0—7 f TRANSITION IN A SINGLE 
JOSEPHSON JUNCTION 


As a first step, let us recall some known results con¬ 
cerning a single dot embedded in a Josephson junction. 
In this section, we summarize the properties of such a 
setup, and we determine under what condition the junc¬ 
tion is in the the 7 r or the O^ 2 ) phase. The 0 
phase is characterized by a positive Josephson current 
for (f) £ [0, 7 r]. It can be further separated between a 
phase where the dot is almost empty and a 0 ^ 2 ^ phase 
where the mean occupation number of the dot is almost 
2. The 7 T phase is associated with a negative Joseph¬ 
son current for (j> £ [ 0 , w], and corresponds to a singly 
occupied dot. 

In perturbative calculations, the Josephson current 
flowing through the quantum dot is the result of the 
tunneling of Cooper pairs, which requires a fourth or¬ 
der perturbative expansion in the tunneling amplitudes 
between a superconductor and the quantum dot. At this 
order, the current can be written as J = Josin^, where 
<f> is the phase difference between the superconductors, 
and the sign of Jq determines the 0 or 7 r character of the 
junction. Providing a single occupancj^n the quantum 
dot (0 < —£ <C U ) Jo can be negative 32 unlike the case 
of an empty quantum dot. 

This phenomenon has been investigated numerically 33 
for arbitrary transmissions between the dot and the su¬ 
perconducting leads. For a fixed quantum dot energy 
level £ < 0 and a fixed decay rate T, the current as a 
function of the phase difference (j) between the super¬ 
conductors undergoes a discontinuity as one tunes the 
Coulomb interaction U across a critical value. Comput¬ 
ing the mean occupation number on the quantum dot in 
the (—£, U) plane reveals the presence of all three phases, 
as shown in Fig. [3] where the 7 r phase, which lies around 
the line 2 e + U = 0, separates the O ^ 0 - 1 and O*' 2 ^ phases. 
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V. SPLITTING EFFICIENCY 


In order to identify the optimal regime for observing 
signatures of CAR processes in the AB signal, we need 
to define a splitting efficiency. We start by discussing the 
case of low electron transmission, where intuition can be 
gained from simple perturbation theory. We then aim at 
comparing the results of our general approach for arbi¬ 
trary transmission to those perturbative results for the 
AB signal. 

The fourth order perturbation expansiorPO - 30 in the 
tunneling amplitudes t :) allows us to write the Josephson 
current as three contributions, associated with the three 
different processes illustrated in Fig. [lj so that 


J(<f>, a) = Id sin (<j> + a) + Iu sin (cj> — a) + Icar sin<^>. 

(23) 

Indeed, in the presence of a magnetic flux a / 0, an ad¬ 
ditional phase shift is acquired, depending on the path 
that each electron of the Cooper pair followed between 
the two superconducting leads. When the Cooper pair 
tunnels through the U (resp. D ) quantum dot, it ac¬ 
cumulates a phase shift —a (resp. +a) in addition to 
the superconducting phase difference and contributes to 
the Josephson current with an amplitude Ijj (resp. Id)- 
However, when the Cooper pair is delocalized on the two 
quantum dots (via the CAR process, thus contributing 
with an amplitude Icar), one electron gets a phase shift 
+ " while the other one gets a phase shift — “ i the pair 
accumulating as a result no additional phase shift. 

The critical current associated with the Josephson cur¬ 


rent Eq. (23) then reads 


J c (a) = I 0 Vl + a cos a + b cos 2 a 


(24) 


with 

and 


P 0 b = 


= Vi? 


u ^ 
2IdIu 


V + ^car> V a = 2/car (Iu + Id) 
For low enough values of T, we are 


in the tunneling regime and the approximation Eq. (23) 
for the Josephson current is justified. If we are able to 
extract the parameters Ijj, Id and Icar, e.g. from a fit 
of the AB signal, we can calculate the quantity 


r t = 


T 2 

-'car 


I? + Id + Var 


(25) 


which encodes the splitting efficiency of the double 
Josephson junction. It varies between 0 and 1. ry = 0 
corresponds to a low efficiency of the CAR process while 
ry = 1 is obtained when this process of spatial delocal¬ 
ization of the two electrons of a Cooper pair is much 
more important than the tunneling processes of a whole 
Cooper pair through a single quantum dot. 

As already stressed out, the formalism developed in 
Secs. m is valid regardless of the strength of the 
coupling between the superconductors and the quantum 
dots. However, the above definition of the splitting ef¬ 
ficiency relies on the expression Eq. (24) for the critical 


current which is no longer valid in the non-perturbative 
regime. There, one needs an alternative diagnosis for the 


detection of the CAR process. As it turns out, a relevant 
quantity can be extracted from the mean powers of the 
critical current obtained for infinitely distant injection 
points (77 = 0) and for coinciding ones (77 = 1). Indeed, 
defining V v = f a * da [J% (a)] 2 , we compute the quantity 


\I~>r)=l — W T7 =01 

Vrj-l 


(26) 


This generalizes the concept of splitting efficiency to the 
case of arbitrary transmission, and coincides with the 
definition of Eq. (25) in the tunneling regime. 


VI. NANOSQUID IN THE TUNNELING 
REGIME 

We first focus on the tunneling regime, taking a low 
value of the decay rate T = 0.01A. We choose to explore 
all the possible combinations for the phases of the two 
quantum dots (see Figs. |4|5]). The results were obtained 
for symmetric couplings 7 = 0, at temperature /I -1 = 
0.002A. The energies of the quantum dots are: £ = 
—0.3A for the 7 r phase, £ = 0.3A for the O ^ 0 - 1 phase and 
£ = —0.9A for the 0^ phase. The Coulomb interaction 
U is chosen to be the same for the two quantum dots, 
and within a specific range as staying in a given phase at 
fixed energy restricts the possible values for U. 

The particle-hole symmetry ensures that the current 
is invariant under the change (eu,£d) —► (—£[/,—£d). 
For the values of energy £ mentioned above, this implies 
that the critical current is identical for 0 ^ — 0 ^ and 
0 (2) — phase associations with U = 0 . 6 A, for tt — 0^ 
and 7T — 0*- 2 ) phase associations again with U = 0 . 6 A, 
and finally for C)( 0 ) — O® at U = 0.7A and — O^ 2 - 1 at 
U = 0.5A 

We first consider the critical current curves of a SQUID 
made of two independent single Josephson junctions, i.e. 



Figure 4. Critical current (in units of 10 4 eA/ft) curves for 
symmetric associations of dots in the tunneling regime: 7r — 7r 
(top panel), 0 (0 ^ — 0 (o) (middle panel), 0^ — O^ 2 - 1 (bottom 
panel). 
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h = °l 



U = 0.5A - U — 0.7A 

V = 0.6A U = 0.8A 


U = 0.9A 
U = l.OA 


Figure 5. Critical current (in units of 10 4 eA//i) curves for 
asymmetric associations of dots in the tunneling regime: 7 r — 
0 <°> (top panel), iv — 0 ^ (middle panel), 0 (0 ^ — O ^ 2 - 1 (bottom 
panel). 


for r] = 0. In this particular case, /car = 0 and con¬ 
sequently a = 0 in Eq. (24), so that the extrema of 
the critical current are the zeros of sin 2a, and the total 
Joseplrson current is given by Eq. ( [23| with /car = 0. 
From the results of Figs. [I] and [5j it clearly appears that 
there is a 7r/2 phase shift in the critical current between 
the situation where the two dots are in the same phase 
(0 — 0 or 7r — 7r) and the one where they are in differ¬ 
ent phases (ir — 0). More specifically, for two quantum 
dots in the same phase ( IdIu > 0), there is no phase 
shift between the currents of the two single Joseplrson 
junctions for a = 0, so that the maxima are added 
(| Id + lu | = \Id\ + |/f/1) and, as a result, the critical 
current is maximal for a = 0. However, for two quantum 
dots in different phases ( IdIu < 0), the phase shift of 
7 r between the currents of the two single Josephson junc¬ 


tions for a = 0 is compensated by a phase shift a = + 7 t /2 
for one of the currents and a phase shift —a = — 7t/2 for 
the other one {\Id ~ Iu\ = \Id\ + \Iu\) resulting in a 
maximum of the critical current for a = 7t/ 2. Such a 
behavior has been observed experimentally.^ 

Comparing the left panels of Figs. [4][5]to the right ones 
(i.e. the case 77 = 0 to 77 = 1 ), we immediately obtain 
evidence of the cross-talk between the two single Joseph¬ 
son junctions: the period of the critical current doubles. 
This is a signature of the emergence of the CAR process. 
Note that the symmetric associations (cf. Fig. [4]) differ 
completely between 77 = 0 and 77 = 1 : for the 7 r — 7 r as¬ 
sociation, the maximum at a = 0 for 77 = 0 becomes a 
minimum for 77 = 1, and for the 0^°^ — 0^°^ and C)( 2 ) — O *- 2 - 1 
associations, the maximum at a = 7 r for 77 = 0 becomes a 
zero for 77 = 1. Concerning the asymmetric associations 
(cf. Fig. [S]), the critical current curves for 77 = 1 somehow 
look like those obtained for 77 = 0 : the positions of the 
maxima and minima are mostly preserved for all phase 
associations, only their local or global character changes 
when tuning 77 . 

Increasing U results in a more pronounced filtering as 
the processes where the quantum dots are doubly occu¬ 
pied are less favored. This explains the observed decrease 
in critical current for the O^ 0 - 1 — 0^ and tt — tt associa¬ 
tions (Fig. [4]). The opposite behavior happens when the 
quantum dots are doubly occupied, i.e. for the 0^ 2 ^ — 0 (2 ^ 
association. There, increasing U favors processes where 
the occupation of the dots is lowered, since approaching 
the 7 T transition results in the decrease of the mean occu¬ 
pation number on the quantum dot. Such opposite be¬ 
haviors while tuning U for 0^°^ —0^°^ and O^ 2 ) — O^ 2 ) phase 
associations can be seen as a consequence of the already 
discussed particle-hole symmetry. Similarly, the 7r — 0^ 
and 7T — O^ 2 - 1 combinations have opposite interaction de¬ 
pendence (see Fig. [ 5 ]). Interestingly, the interaction has 
no noticeable effect in the case of the 0^°^ — 0^ phase 
association. 
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Figure 6 . Influence of the parameter 77 on the critical current 
(in units of 10 4 eA/?i) in the tunneling regime (U = 0.6A). 


Figure 7. Splitting efficiency r given by Eq. (26 1 extracted 
from the curves of Figs. [3] and [5] 





























In order to monitor the emergence of CAR processes, 
we investigate in Fig. [^intermediate regimes for the tt—tt, 
0 ( 2 ) — o^ 2 ), 7 r — 0< 2 > and — C/ 2 ) phase associations 
by turning on 77 progressively from 0 to 1. This is a 
phenomenological way to introduce more and more CAR 
processes, as the injection points in the superconductors 
are brought together, from 77 = 0 (no CAR effect) to 
77 = 1 (maximal CAR effect). We observe a dramatic 
increase of the critical current in the symmetric configu¬ 
rations (tt — n and O*- 2 ) — C)( 2 ) phase associations), while 
in the asymmetric configurations (n — O ^ 2 - 1 and C/ 0 * — C/ 2 ) 
phase associations), the modification of the AB signal is 
noticeable, but not substantial. 


We display the splitting efficiency r given by Eq. (26) 
for all phase associations in Fig. [7} In the studied domain 
for Coulomb repulsion parameter U, we do not observe 
noticeable evolutions of r except for the — n associa¬ 
tion for which we see a clear decrease. The symmetric as¬ 
sociations of phases (for which the critical current differs 
the most between 77 = 0 and 77 = 1 ) lead to the highest 
values of r. The highest splitting efficiency is obtained 
for two quantum dots in the 7 r phase: ensuring a mean 
occupation number around 1 on each quantum dot favors 
the CAR process. One can get some insight concerning 
the value of r for the associations O® — O ^ 0 - 1 and O^ 2 -* — C/ 2 ) 
from the perturbative calculation presented in the previ¬ 
ous section. Indeed, in the tunneling regime the splitting 
efficiency r matches the definition r t of Eq. (25) in terms 


of the amplitudes Id, lu and /car- For symmetric asso¬ 
ciations, we must add the constraint Id = Ijj = I /2 to 
the fit procedure, so that the critical current reads 


J c {a) 

1*1 


/car 


+ COSO! 


(27) 


and since the critical current in symmetric associations 
of 0 phases vanishes for a = n, we get /car// = 1 and 
r t = 2/3. 


to reproduce the qualitative study of Sec. [VI] Our goal is 
to determine what features are preserved and what has 
changed, and to compute the splitting efficiency r defined 
by Eq. (26) in order to determine which associations of 


phases favor nonlocal phenomena the most. The critical 
current curves are given in Figs.[8][9j 

From the results of the tt — tt association for 77 = 1, it is 
clear that we are no longer in the tunneling regime as the 
flux dependence can not be fitted by Eq. (24) together 


with the constraint that Id = lu■ We can again notice, 
for 77 = 0 , the 7t/2 phase shift in the critical current 
between dots in the 7r — 0 phases (top and middle panels 
of Fig. [9]) and dots in the 0 — 0 or 7r — 7r phases (Fig. [8] 
and bottom panel of Fig. [9]). 

The critical current in Figs. [§][ 9 ] again presents a dou¬ 
bling of its period when switching 77 from 0 to 1 . This is 
due to the emergence of nonlocal processes where both 
quantum dots are involved. For the 7r — 7r association, 
during this switching, a = 0 remains a local maximum 
(contrary to the tunneling regime). For O® — C/ 0 ) and 
(7 2) — C/ 2 ) phase associations, the maximum at a = tt 
for ?/ = 0 becomes a zero for 77 = 1 (as in the tunneling 
regime). There is little influence (less than in the tunnel¬ 
ing regime) of 77 on the 7r — and n — C/°) associations 
whereas, for the O*- 0 ) — Ch 2 ' association the maximum in 
a = 0 is considerably lowered (more than in the tunneling 
regime) from 77 = 0 to 77 = 1. 

For symmetric associations of phases, the evolution of 
the critical current with increasing U (Fig. [sj) can be ex¬ 
plained following the same arguments as in the tunneling 
regime. While the filtering of electrons tunneling through 
0 (o) or tt quantum dots is responsible for a decrease of 
the signal, favoring one-electron processes through a 0*- 2 -* 
quantum dot results in an increase of the signal. The 
opposite 17 -dependence for the 0*- 0 -* — C)( 0 ) and C)( 2 ) — 0*- 2 ) 
phase associations is a consequence of particle-hole sym- 


VII. NANOSQUID IN THE HIGH 
TRANSPARENCY REGIME 

The advantage of the formalism developed in Secs. |TT[ 
|III| relies on the possibility to address high transparency 
regimes of the double Josephson junction under study. 
The results presented in this Sec. are obtained for sym¬ 
metric couplings 7 = 0, at temperature /3 _1 = 0.02A 
and for a decay rate T = 2A. The energies of the quan¬ 
tum dots are: e = —4A for the tt phase, e = 4A for the 
0(°) phase and e = —12A for the f/ 2 ) phase. Particle- 
hole symmetry implies that the critical current is identi¬ 
cal for 0 (°) — 0 <°> and C/ 2 ) — C)( 2 ) phase associations with 
U = 8 A, for 7 T — 0*- 0 ) and 7 r — 0 ( 2 ^ phase associations again 
with U = 8A, and finally for 0 (o ) - C>( 0 ) at U = 9A and 
0 ( 2 ) - o( 2 ) at U = 7A. 

The strategy is again to investigate all the possible 
combinations of phases for the quantum dots (which cor¬ 
respond to different mean occupation numbers), in order 




0 0.5 1 1.5 2 0 0.5 1 1.5 2 

I 771 <2/7r a /tt \ 71 

77 = 0 ' 1 \rj = 1 

U = 7A - U = 9A 

U = 8A U = 10A 

Figure 8. Critical current (in units of eA/h) curves for sym¬ 
metric associations of dots in the high transparency regime: 
tt—tt (top panel), 0 (O ^—O^ 0 - 1 (middle panel), 0 ( A— (/ 2 ) (bottom 
panel). 



























9 



U = 7A - U = 9A 

U = 8A . U = 10A 


Figure 9. Critical current (in units of eA/h) curves for asym¬ 
metric associations of dots in the high transparency regime: 
7T—0^°' (top panel), n — O ^ 2 - 1 (middle panel), 0 lH) —O ^ 2 - 1 (bottom 
panel). 


metry. While there is still no noticeable effect of U on 
the association as in the tunneling regime, the 

7 r — association also shows little 17-dependence. As a 

consequence, we do not observe the opposite behavior as 
a function of U for the 7 r — (A 0 ^ and 7 r — 0 (2 ^ associations 

(Fig. [o]). 

We introduce progressively nonlocal effects in Fig. [To] 
by switching i] from 0 to 1. For the symmetric associa¬ 
tions ( 7 T — 7 r and — 0 i°i) as well as for the 0 ^ — 0 *- 2 ) 
configuration (contrary to the tunneling regime), the AB 
signal is dramatically increased when CAR processes are 
switched on. On the contrary, the 7 r — 0^ 2 i phase asso¬ 
ciation is hardly influenced by the presence of nonlocal 
processes (even less than in the tunneling regime). We 
can quantify the importance of CAR processes compared 
to the direct tunneling through a single quantum dot by 


calculating the splitting efficiency r given by Eq. (26) 


which we display in Fig. |TT| for the different associa¬ 
tions of phases. As a function of U, the splitting effi¬ 
ciencies which are found are essentially constant over the 
range considered. As it turns out, although we work with 
specific phases for the individual quantum dots, imply¬ 
ing specific populations, the quantization of the electron 
charge on these dots is ineffective at high transparencies 
because they constitute “open quantum dots”, involving 
large fluctuations from their average population. This 
is consistent with the observation that the splitting effi¬ 
ciency varies little with U. As mentioned above, there is 
little influence of 77 on the 7r— O^ 2 ) and 7r— O*- 0 ) associations 
and this is why these are the lowest splitting efficiencies 
we found. The splitting efficiencies of the associations of 
0 phases are around 0.55 and the highest splitting effi¬ 
ciency is still obtained for the 7 r — 7 r association. 

Finally, we note that the energy levels of the dots can 
be easily varied experimentally using electrostatic gates. 
Varying the position of the energy level £ of a quantum 
dot at fixed decay rate T and fixed Coulomb on-site re¬ 
pulsion U changes the effective transparency defined as 


D(e) = 


r 2 +( £ +^r 


(28) 


Thus, by varying the energy levels of the two dots, we 
can optimize the splitting efficiency by reaching the n — 7 r 
phase at high effective transparency. We show on the top 
panel of Fig.[l2]the splitting efficiency as a function of the 
energy level of the two dots, taken to be identical. The 
abrupt change of the efficiency near D(e) = 0.4 shows 
the crossover from the 0 ^°^ — 0 ^°^ to the n — n junction. 
Similarly, when tuning only one of the two energy levels, 
while maintaining the other dot in the 7 r (middle panel) 
or the 0 ^°) phase (bottom panel), the splitting efficiency 
shows a marked transition as a function of the effective 
transparency. 
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Figure 10. Influence of the parameter 77 on the critical current 
(in units of eA/h) in the high transparency regime (U = 8 A). 


Figure 11. Splitting efficiency r given by Eq. (26 1 extracted 
from the curves of Figs. [ 8 ] and [9] 
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4 3 2 1 Q -1 -2 -3 -4 



Figure 12. Evolution of the splitting efficiency as a function 
of the effective transparency D(e) [Eq. ( |28| )], for fixed F = 
2A and U = 8A. Top: The energy levels of both dots are 
taken to be identical and varied simultaneously ( eu = £d = 
e). Middle: the D quantum dot is taken in the 7r phase 
(sd = —4A) while the energy of the U dot is varied (eu = e). 
Bottom: the D quantum dot is taken in the 0 (o) phase (eo = 
4A) while the energy of the U dot is varied (eu = e). 


VIII. CONCLUSION 

We have studied a double Josephson junction consist¬ 
ing of two quantum dots connected to two superconduc¬ 
tors as a tool to probe Cooper pair splitting. When the 
two single Josephson junctions which constitute our de¬ 
vice are coupled to each other via CAR processes, the 
doubling of the period of the critical current measured 
in a nanoSQUID experiment is an evidence for the emer¬ 
gence of nonlocal phenomena in the electronic transport 
of the double junction. This type of diagnosis may prove 
more convenient than non-equilibrium scenarios involv¬ 
ing a superconducting source of electrons and two nor¬ 
mal leads/dots where either nonlocal conductance signal 
or noise cross-correlations are measured. 

While this device had been studied in the context of 
perturbation theory in the tunneling Hamiltonian cou¬ 
pling the dot to the leads 1321321 n0 generalization to ar¬ 
bitrary transmission had been proposed, and no system¬ 
atic study for optimizing the CAR signal with respect 
to the phases (O®, n, CK 2 )) of the quantum dots had 
been attempted so far. The path integral approach of 
Ref. l33l within reasonable approximations (saddle point 
treatment) allows precisely to meet these goals. One of 
the key results of the present work resides in defining the 
degree of efficiency of Cooper pair splitting, and evaluat¬ 
ing it for the different dot configurations. 

We first studied the tunneling regime where the usual 
perturbative expression for the Josephson signal allows 
us to fit our numerical results and to introduce a nat¬ 
ural definition of the splitting efficiency. We were also 
able to generalize this quantity to arbitrary transparency, 


providing a criterion for the efficiency of CAR processes 
which is based on an analysis of the AB signal of the 
Josephson critical current. We thus studied the promi¬ 
nence of nonlocal phenomena depending on the phases 
of the quantum dots and found that the 7r — 7r associa¬ 
tion optimizes the splitting of the Cooper pairs that are 
emitted from one superconductor and recombined on the 
other one. Yet, our analysis shows that the 0® — 0^ 
and (b 2 ) — O^ 2 - 1 combinations also provide robust Cooper 
pair splitting signals. Within each of these combinations 
of phases, we see for the most part that variations of the 
Coulomb repulsion parameter has little influence on the 
Cooper pair splitting efficiency. 

The present treatment should only be valid if 
the su perco nductor-dots system is above the Kondo 
regime 1321321 F or a single dot embedded between two su¬ 
perconductors, the Kondo effect manifests itself when 
the Kondo temperature is larger than the superconduct¬ 
ing gap, because quasiparticle excitations are required 
to trigger the spin flip between the leads and the dot. 
The expression of the Kondo tempera ture for a Hub¬ 
bard type Coulomb repulsion is giveiP^S! by t k = 
VTU /2 e ne o(so+u)/(ru) _ Within the range of parameters 
chosen in our numerical study, we find that the Kondo 
temperature is at most 0.13A, which gives us confidence 
in our working assumptions. At any rate, the Kondo 
regime could be avoided by working with the 0*°' — 0 (o ) 
combination of phases, which according to Fig. |TT] still 
has a sizable splitting efficiency. 

Furthermore, we treated the CAR coupling parameter 
r) in a phenomenological manner (varying it between 0 
and 1), justifiably so because so far no experimental in¬ 
vestigation of Cooper pair splitting in superconducting- 
normal metal “forks” out of equilibrium seems to find the 
power law suppression which is attributed to bulk 3D su¬ 
perconductors. This may be due to the fact that micro¬ 
scopic models have to be revisited, taking into account 
that electron emission/absorption in the superconductors 
is not point like, but should be averaged over some finite 
volume, reducing the effect of oscillations over the Fermi 
wavelength of the CAR parameter. 

Extensions of this work could be envisioned by going 
beyond the saddle point approximation of the Hubbard 
Stratonovich transformation. This transformation is ex¬ 
act in its functional integral form and neglecting the fluc¬ 
tuations of the auxiliary field and then using a method of 
steepest descent are approximations which were sufficient 
to exhibit the 7r-shift in a single junction^. However, the 
possibility to go beyond this mean field type of approach 
could be considered by looking at Gaussian fluctuations 
around the stationary value of the auxiliary held. Al¬ 
ternatively, a numerical renormalization group method 38 
could be employed in principle, but the proliferation of 
couplings and parameters is likely to render it cumber¬ 
some. 

A limitation of the present Cooper pair diagnosis re¬ 
sides in the fact that the nanoSQUID requires reduced 
dimensions so that the injection points to the dots on 
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Figure 13. Setup with nanowire/nanotube quantum dots. 

both superconductors are separated within a distance 
smaller than the coherence length. This imposes con¬ 
straints on the separation of the quantum dots, and as 
a consequence, the area where the AB flux is imposed 
becomes reduced. Recall that in order to perform the 
AB diagnosis, a few flux quanta need to be introduced 
in this area, but if the imposed magnetic field needed to 
apply several flux quanta becomes larger than the critical 
field of the superconductors, the whole diagnosis breaks 
down. In order to optimize the surface area encompassed 
between the dots, we believe the best choice would be 
to work with nanowire/nanotube quantum dots, which 
have a large aspect ratio, and which nevertheless can 
achieve large charging energies even though their length 
may exceed several /im. 41 This possibility is illustrated 
in Fig. |13[ To work out the numbers, we find that im¬ 


posing 2 magnetic flux quanta within a 1 /rm 2 area re¬ 
quires a magnetic field of 8 x 10~ 3 Tesla, which is still 
smaller than the critical field of superconducting mate¬ 
rials such as Aluminum (10~ 2 Teslas) or Niobium (0.2 
Teslas). Note that Aluminum could be quite adapted 
due to its long coherence length (1.6 /jm). 

While this work came to completion, we became aware 
of a current biased measurement 42 which studies pre¬ 
cisely the behavior of the double dot Josephson junction. 
There, a (remarkable) comparative study of the switch¬ 
ing current (the current required to transit to the dissi¬ 
pative regime with voltage bias) was performed for differ¬ 
ent dot phase configurations. This experiment seems to 
provide realistic evidence of Cooper pair splitting, albeit 
indirectly, because the self organized quantum dots em¬ 
bedded in the Josephson junction are too close in order 
to impose the necessary magnetic flux to observe the AB 
oscillations. 
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